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ABSTRACT 

This is the first paper of a series wfiere we study tlie ciustering of LRG galaxies 
in the latest spectroscopic SDSS data release, DR6, which has 75000 LRG galax- 
ies covering over 1 Gpc? /h^ at 0.15 < 2 < 0.47. Here we focus on modeling red- 
shift space distortions in ^(tt, cr), the 2-point correlation in separate line-of-sight and 
perpendicular directions, on large scales. We use large mock simulations to study 
the validity of models and errors. We show that errors in the data are dominated 
by a shot-noise term that is 40% larger than the Poisson error commonly used. We 
first use the normalized quadrupole for the whole sample (mean z=0.34) to estimate 
13 = f(p.rn)lb = 0.34 ± 0.03, where f{^m) is the linear velocity growth factor and 
b is the linear bias parameter that relates galaxy to matter fluctuations on large 
scales. We next use the full Hit, a) plane to find f^om = 0.245 ± 0.020 (h=0.72) and 
the biased amplitude hag, = 1.56 ± 0.09. For standard gravity, we can combine these 
measurements to break degeneracies and find (j% ~ 0.85 ± 0.06, 6 — 1.85 ± 0.25 and 
f{flm) — 0.64 ± 0.09. We present constraints for modified theories of gravity and find 
that standard gravity is consistent with data as long as 0.80 < as < 0.92. We also 
calculate the cross-correlation with WMAP5 and show how both methods to measure 
the growth history are complementary to constrain non-standard models of gravity. 
Finally, we show results for different redshift slices, including a prominent BAO peak 
in the monopole at different redshifts. The ^(tt, a) data on large scales is shown to 
be in remarkable agreement with predictions and shows a characteristic large region 
of negative correlation in the line of sight, a BAO ring and a prominent radial BAO 
peak. The significance of this is presented in paper IV of this series. We include a 
study of possible systematic effects in our analysis to find that these results are quite 
robust. 



1 INTRODUCTION 



Galaxy clustering allows the study of different physical phe- 
nomena at each scale. On large scales density fluctuations 
are small and can be modeled by linear theory to constrain 
cosmological parameters. We need large surveys with galaxy 
positions to do this. 

The measured redshift distance s of a galaxy differs from 
the true radial distance r by its peculiar velocity along the 
line-of-sight «,-. These displacements lead to redshift dis- 
tortions, with two important contributions. The first, on 
large scale fluctuations, is caused by coherent bulk motion. 
We see walls denser and voids bigger and emptier, with a 
squashing effect in the 2-point correlation function along the 
line-of-sight, known as the Kaiser effect (Kaiser 1987). At 
small scales, random velocities inside clusters and groups of 



galaxies produce a radial stretching pointed at the observer, 
known as fingers of God (FOG) (Jackson 1972). 

Although such distortions complicate the interpretation 
of redshift maps as positional maps, they have the advantage 
of bearing unique information about the dynamics of galax- 
ies. In particular, the amplitude of the quadrupole redshift 
distortion on large scales yields a measure of the linear red- 
shift distortion parameter /3, which is related to the matter 
density, and gives direct information on the growth of the 
newtonian gravitational potential /(fi) and bias b, which 
describes how the number density of galaxies traces matter 
density (see Eq.(|3])). We will also show how one can use the 
anisotropic redshift distortions, broken into line of sight and 
perpendicular directions, to find separate measurements for 
/(fi) and b. 

In this work we use the most recent luminous red galax- 
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ies (LRGs) from the spectroscopic SDSS public data release 
DR6 (Adelman-McCarthy et al 2006), and perform studies 
of linear bias on large scales to obtain cosmological param- 
eters. As we will see, we can break the degeneracy between 
bias and as present in the correlation function thanks to red- 
shift distortions anisotropics , and look at the growth history 
and possible modifications of the gravity (see Linder 2005, 
Guzzo et al 2008 and the recent Yamamoto et al 2008). 

A dark energy dominated universe causes gravitational 
potentials to evolve, changing the energy of CMB photons 
passing through them. This is the so called Integrated Sachs- 
Wolf (ISW) effect. Here, we have also cross-correlated LRGs 
with WMAP (Hinshaw et al 2008) in order to investigate 
this ISW effect, reproducing a high signal as seen in recent 
studies (see Gaztanaga etal 2006, Giannantonio et al 2008), 
as compared to current A.CDM model. In this way, we can 
also break the bias — ag degeneracy and study the growth 
history using the cross- correlation between temperature of 
CMB and galaxies wtg , which give independent constraints 
on bias and ag, compared with those obtained from redshift 
space distortions. 

The same LRGs (but with reduced area) have been 
studied from different points of view. Tegmark et al (2006) 
have done an analysis of the power spectrum at large scales 
to obtain cosmological parameters. Zehavi et al (2005) study 
LRGs at intermediate scales (0.3 to 40Mpc/h), where they 
calculate the projected correlation function, the monopole 
and real-space correlation function to study mainly the lin- 
ear high bias, the non-linear bias and the differences be- 
tween luminosities, remarking that there are differences from 
a power law for scales smaller than IMpc/h. 

While we were working on our results, a very interesting 
paper by Okunmra et al (2008) showed the first application 
of the anisotropy in the 2-point correlation function includ- 
ing the baryonic features, to obtain constraints on cosmo- 
logical parameters. They use the DR3 spectroscopic sample 
of LRGs to calculate the 2-point correlation function with a 
different definition of what parallel and perpendicular direc- 
tions mean (the line of sight is defined in their paper as the 
direction toward the galaxy 1, Matsubara 2004). They fit 
^(<T, tt) for large scales using the linear Kaiser model, from 
40Mpc/h to 200Mpc/h, excluding the FOG zone. They point 
out that a direct measurement of the growth function can be 
obtained from the distortion to the clustering signal caused 
by bulk flows (Kaiser anisotropy). Although constraints are 
weak by now, the fltting of all the anisotropic 2-point corre- 
lation function, including the baryonic feature, would enable 
us to divide the effect of redshift distortions into dynamical 
and geometrical components. The anisotropy due to geo- 
metric distortion contributes to a better estimation of the 
equation of state of the dark energy. Here we try to obtain 
new constraints, encouraged by the fact that DR6 has dou- 
bled the volume of DR3. 

There has already been significant progress in this direc- 



tion using LRG. Eisenstein et al (2005) detected the baryon 
acoustic peak in the 2-point correlation function using LRGs. 
Hiitsi (2006 a,b) use LRGs to constrain cosmological param- 
eters in the power spectrum, including the baryonic peak. 
Percival et al (2007) have analyzed also the LRGs using both 
2dF and SDSS. Padmanabhan et al (2007) used the photo- 
metric catalog to work with a bigger set of LRGs with pho- 
tometric redshifts, obtaining also cosmological constraints, 
the same as Blake et al (2007), which work with the MegaZ- 
LRG, a photometric-redshift catalog of luminous red galax- 
ies based on the imaging data of the SDSS DR4. 

This is the first and main paper of a series on clustering 
of LRG. Here we will focus on redshift space distortions in 
^(7r,a) on larges scales. In Paper II (Cabre and Gaztanaga 
2008) we will look at small scales in ^(tt, a) and address the 
issue of the FOG. In Paper HI (Gaztanaga et al 2008) we 
look at the 3-point function and in Paper IV (Gaztanaga 
et al 2008) we will focus on the significance of the BAO 
detections. In this first paper we present the basis for the 
error analysis and the study of systematics effects that will 
be used in the rest of the series. 

The paper is organized as follows. Section II gives a 
summary of the theory involved in redshift space distor- 
tions, the model used for the 2-point correlation function 
in separate line of sight and perpendicular directions, and 
the different multipoles. Section HI presents an accurate de- 
scription of the data selected to work with, LRG galaxies. 
Section IV is the main part of this paper, where we explain 
the analysis done and our main results. We use the distor- 
tions in redshift space to obtain the distortion parameter 13, 
Q,m, OS and modifications to gravity. We also look at the 
ISW effects for this LRG data, and explore the differences 
between redshift slices. As usual, we end with discussion 
and future work in section V. The Appendix A is devoted 
to present the mock simulations that we use to validate the 
errors used, Monte Carlo, jackknife and a new theoretical ap- 
proach. In this section, we also validate the methods used to 
obtain constraints (section IV) in real data with our mocks. 
In the Appendix B we present a study of possible systematic 
effects. 



2 THEORY 

Kaiser (1987) pointed out that, in the large-scale linear regime, 
and in the plane-parallel approximation (where galaxies are 
taken to be sufficiently far away from the observer that the 
displacements induced by peculiar velocities are effectively 
parallel), the distortion caused by coherent infall velocities 
takes a particularly simple form in Fourier space: 



Ps{k) 



{l + PdfP{k). 



(1) 



where P{k) is the power spectrum of density fluctuations 5, 
H is the cosine of the angle between k and the line-of-sight. 
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the subscript s indicates redshift space, and /3 is proportional 
to the velocity growth rate in linear theory. 

If the galaxy overdensity 5 is linearly biased by a factor b 
relative to the underlying matter density 5m of the Universe, 

S = b5m, (2) 

then the observed value of l3 is 

fiflm) _ IdlnD 



and 



P- 



(3) 



b b din a 

where D is the linear density growth factor and / the veloc- 
ity growth rate, which we can write as: 

f{Q.m) = ^m{ay (4) 

where 7 is the gravitational growth index (Linder 2005), 
Q,Tn{a) is the matter density at a redshift z where a = 1/(1 + 



^m{a) — 






(5) 



where fiom is the matter density at z=0 (which is often 
just called £!„) and 



H{a) = Ho ^noma-'^ + (1 - nom)a-3(i+'") (6) 

where w is dark energy equation of state parameter. 

By construction (Linder 2005), the growth index for- 
malism separates out two physical effects on the growth of 
structure: fl{a) involves the expansion history and 7 focuses 
on the gravity theory. The value 7 = 0.55 corresponds to 
standard gravity, while 7 is different for modified gravity, 
for example 7 = 0.68 in the braneworld cosmology (Linder 
2005). The linear density growth factor can be found using 
Eq.©: 



D(a) = e 



f^ dlna [n{a)"'-ll 



(7) 



2.1 Modelling ^(CT,7r) 

Hamilton (1992) translated Kaiser's results into real space, 

C'(a,7r) =eo(s)Po(M)+6(s)^2(/i)+^4(s)P4(M), (8) 
where tt is the separation along the line-of-sight (LOS) and 
a is the separation in the plane of the sky ,the absolute 
distance of separation is s = y/a'^ + n^, fi is the cosine of the 
angle between s and the line-of-sight and Pt are Legendre 
polynomials. 

In general, we can define the multipoles of ^(vr, a) (see 
also Eg. fTell ). 



6(^)=(Y + ^)K(r)-C(r)], 



Us) 



35 



C(0 + l^r) - l^r) 



(9) 

(10) 

(11) 



^(^) = -T / £.{r'y^dr', 



S.ir) = —, I £.{r')r'*dr'. 



(12) 



(13) 



We use these relations to create a model 5' (a, n). In this 
paper, which analyzes large scales, we will use a linear real 
space correlation function ^(r) as input (Eisenstein and Hu 
1998) u Then we convolve it with the distribution function 
of random pairwise velocities, f{v), to give the final model 
C(cr,7r) (Peebles 1980): 



a<^, ^) = / C'(^, ^ - v/Hiz)/a{z))fiv)dv 



(14) 



where we divide peculiar velocities by a{z) to translate 
to comoving distances, since velocities are defined in physical 
coordinates. 

We represent the random motions by an exponential 
form, 



/(«) 



(T„\/2 



exp 



V2U 



(15) 



where (t„ is the pairwise peculiar velocity dispersion. An 
exponential form for the random motions has been found 
to fit the observed data better than other functional forms 
(Ratcliffe et al 1998; Landy 2002). We have also checked 
this relation with our simulations, which will be described 
later on. 

Matsubara (2004) and Scoccimarro (2004) have pre- 
sented different models for the 2-point correlation function 
in redshift space. Tinker et al (2006) and Tinker (2007) 
model redshift space distortions in the context of halo oc- 
cupation distribution model (HOD). These models are com- 
plementary to the one studied here. In most situations the 
differences are small and we will show that our modeling 
gives good agreement with simulations and real data. 

2.2 Multipoles of ^(a, tt) 

We can define the multipoles of ^{TT,a) as 

6(s) = ^— y a^,a)Pdf^)dfi. (16) 

where /i is cosine of the angle to the line-of-sight n. 
The normalized quadrupole is defined as (Hamilton 1992) 



Q(s) = 



Us) 



(17) 



^o{s)-{3/s^)J^^o{s')s'^ds' 

Linear bias in Eq.([2)) cancels in the quadrupole, and 
Q{s) is only very slightly dependent on the shape of the 



^ We have also tested the inclusion of the effect of non-linear 
growth (with the halofit model) and non-linear bias. We find no 
changes on our fit to large scales, as expected. 
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correlation function, such as changes on large scales due 
to varying Q.m or scale dependent non-linear bias for small 
scales. 

In the Kaiser approximation, ie at large scales, the quadrupole 
is directly related to 13: 



Q{s)^ 



1 + 



(18) 



1/3 +1/32 

The validity of this approximation is studied using simula- 
tions in section |X] At small scales, the quadrupole depends 
strongly on the random pairwise velocities, represented by 
o"„, but it does not depend much on non- linear bias, as we 
will show. 

2.3 The real-space correlation function ^(r) 

We can estimate the real-space correlation function by calcu- 
lating the projected correlation function, H((t), integrating 
the redshift distorted ^(ct, it) along the line-of-sight tt. 



H(a) = 2 



^((T, 7r) dyr 



(19) 



We would like iimax ~ oo, however, with real data, we 
can not integrate until infinity. Here we will use TTmax = 
80Mpc/h. The result does not change when we change the 
upper limit of the integral for TTmax > 60AIpc/h and large a 
in the data. 

Davis and Peebles (1983) show that H((j) is directly 
related to the real-space correlation function. 



H(.) = 2 r ^iM^ 

J IT (r'-^ — a^) 2 



(20) 



It is possible to estimate ^(r) by directly inverting :z{a 
(Saunders et al 1992) 

1 r°° (dE{a)/da 
n 



ar) 



-da. 



(21) 



(a2-r-2)3 

Assuming a step function for H(a) = Hi in bins centered on 
ai, and interpolating between values. 



1 



for r — ai 



^ 7>1 ""'+' ^ ""' 



=i In I ^, I (22) 



CT, + 



y^ 



2-cr2 



We take the redshift space anisotropic model ^(cr, tt) in 
Ea. (|14p . with a fixed j3 — 0.35 and a^ = 400km/ s, and we 
use Ea. (|19|) to obtain the projected correlation function H(cr) 
as if it was data. We use different values to fix the upper limit 
in the integral: n^ax = 60Mpc/h, 80Mpc/h, lOOMpc/h and 
200Mpc/h (color lines in top panel of Fig[T]), and compare 
the "true" result obtained from Eq. (|20|) (solid line in top 
panel of Fig[l}. As we increase n,nax we approach the true 
result, but we can not integrate until 200Mpc/h in real data 
because, as will be shown in the following sections, data 
is quite noisy for n > 60Mpc/h. The implication of these 



100.00 r^ 



10.00 r 



-1.00 r 



0.10 r 



20 40 60 80 100 120 

r(Mpc/h) 




Figure 1. Top: Solid line shows the projected perpendicular 
correlation function S((t) calculated theoretically from ^(r) in 
Eq.lllOj, for a model with /3 = 0.35 and o-„ = 400fcm/s. This 
is compared to integral of ^(tt, cr) in Ea. 11191 1 with 7rmaa;=60.(red 
dashed), 80. (blue dashed), 100. (red dash-dot), 200. (blue dash- 
dot). Bottom; Estimation of the real space correlation function by 
deprojecting the color lines of S(r) in the top panel using Ea.l l21ll . 
We can see good agreement on scales smaller than 30Mpc/h. 



results is that we will not be able to recover the true values 
of S,{a) from S^{-k, a). But we can use this method at smaller 
scales (see Paper II). 

In the bottom panel of Fig[T]we see the real-space cor- 
relation function recovered from the previously calculated 
projected correlation function with different nmax, Eq. (|2ip . 
We obtain a good estimation of ^(r) below 30Mpc/h, where 
we will study the non-linear bias (see Paper II) . 

This analysis has been done for different values of /3 and 
av and we find very similar conclusions. For illustration, we 
have shown in our plots the values that are more in concor- 
dance with real LRG SDSS data. 

Once we recover the real-space correlation function, we 
can also estimate the ratio of the redshift-space correlation 
function, ^(s), to the real-space correlation function, S,{r), 
which gives an estimate of the redshift distortion parameter, 
/?, on large scales: 



C(r) -^ 3 ^ 5 • 



(23) 
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Figure 2. Comoving density vs redshift for our selected LRGs 

As we have shown in Fig.((T}, ^(r) will be in general 
slightly overestimated at large scales when we estimate it 
from the projected correlation function H(cr), so the expres- 
sion jA will in general be slightly lower than expected on 



large scales. 



3 THE DATA 

The luminous red galaxies (LRGs) are selected by color 
and magnitude to obtain intrinsically red galaxies in Sloan 
Digital Sky Survey (SDSS). See Eisenstein et al (2001) or 
http://www.sdss.org for a complete description of the color 
cuts. In this paper, we work with the last public catalog, 
DR6, which covers a solid angle of 47r/7 in the sky. These 
galaxies trace a big volume, around IGpc^ h~^ , which make 
them perfect to study large scale clustering. LRGs are sup- 
posed to be red old elliptical galaxies, which are usually 
passive galaxies, with relatively low star formation rate. 

They have steeper slopes in the correlation function 
than the rest of galaxies, since in some cases, there is more 
than one LRG per halo. They trace a big volume, so they 
represent good candidates to use as dark matter clustering 
tracers. 

LRG's are targeted in the photometric catalog, via cuts 
in the (g-r, r-i, r) color-color-magnitude cube. Note that all 
colors are measured using model magnitudes, and all quan- 
tities are corrected for Galactic extinction following Schlegel 
et al (1998). The galaxy model colors are rotated first to a 
basis that is aligned with the galaxy locus in the (g-r, r-i) 
plane according to: 

cx= (r-i) - (g-r)/4 - 0.18 

C||= 0.7(g-r) + 1.2[(r-i) - 0.18] 

Because the 4000 Angstrom break moves from the g 
band to the r band at a redshift z ~ 0.4, two separate sets 
of selection criteria are needed to target LRGs below and 



above that redshift: 

Cut I for z < 0.4 

rPetro < 13.1 -h C|| / 0.3 

rPetro < 19.2 

|cxl < 0.2 

mu5o < 24.2 mag arcsec~^ 

rpSF - Tmodel > 0.3 

Cut II for z > 0.4 

rPetro < 19.5 

|cxl > 0.45 - (g-r)/6 
g-r > 1.30 + 0.25(r-i) 
muso < 24.2 mag arcsec~^ 

rpSF - Tmodel > 0.5 

where muso is the mean surface brightness within the 
radii containing 50% of Petrosian flux. 

Cut I selection results in an approximately volume- limited 
LRG sample to z=0.38, with additional galaxies to z ~ 0.45. 
Cut II selection adds yet more luminous red galaxies to z ~ 
0.55. The two cuts together result in about 12 LRG targets 
per deg^ that are not already in the main galaxy sample 
(about 10 in Cut I, 2 in Cut II). The radial distribution 
and magnitude- redshift diagrams are plotted in Fig lBll and 

Figim] 

We k-correct the r magnitude using the Blanton pro- 
gram 'kcorrect' Q We need to k-correct the magnitudes in 
order to obtain the absolute magnitudes and eliminate the 
brightest and dimmest galaxies. We have seen that the pre- 
vious cuts limit the intrinsic luminosity to a range —23.2 < 
Mr < —21.2, and we only eliminate from the catalog some 
few galaxies that lay out of the limits. Once we have elimi- 
nated these extreme galaxies, we still do not have a volume 
limited sample at high redshift but we will account for this 
using a random catalog with identical selection function. 

We show the comoving density in Fig[2] once we have 
removed the brightest and dimmest galaxies, and in Fig[3] 
we show the redshift space distribution in a slice of dec = 
32-40 deg. We compare it with the same slice of random 
points. We see clear evidence of clustering in the data. 

We have masked the catalog using at the flrst step the 
photometric DR6 mask, based on the number of galaxies per 
pixel. In previous works (Cabre et al 2006) we showed that 
the mask that we obtain statistically by dropping out the 
pixels with small number of galaxies (see §B3 for details) 
gives identical correlation function that the one obtained 
when extracting the polygons masked by the SDSS team. 
After that, we compare our masked catalog to the LRG 
spectroscopic catalog, and we extract the galaxies that lay 
outside from the "good" plates. Fig |B9l illustrates the result. 

The mask could imprint spurious effects at very small 



http://cosmo.nyu.edu/blanton/kcorrect/kcorrect_help.html 
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DATA 




RANDOM 



z=0.1 



z-0.4 



z=O.S 




Figure 3. Slice in dec = 32-40 deg showing ra vs redshift, in order to see the structure of the data (left panel). We also plot the random 
distribution for this slice (right panel) 



scales. But we are not focused in such small scales, because 
fiber collisions in the redshift catalog are limiting our anal- 
ysis. This is for distances less than 55arc sec, ie less than 
0.3Mpc/h at the mean redshift of LRG data, z=0.34 (see 
Fig. 17 below). We obtain 75,000 galaxies for the final cat- 
alog, from z=0.15 to z=0.47, ie over one cubic Gigaparsec 



4 RESULTS 

In this section we analyze the LRGs from SDSS DR6. First, 
we have used credible simulations to study the errors care- 
fully as well as the methods that we use here (see Appendix 
A for an extensive description). 

Here we define the parameters that we assume during all 
this work, which are motivated by recent results of WMAP, 
SNIa and previous LSS analysis B n^ = 0.98, Q,b = 0.045, 
h = 0.72. We will use the power spectrum analytical form 
for dark matter by Eisenstein and Hu (1998), and the non- 
linear fit to halo theory by Smith et al (2003). For part of our 
analysis we have followed the method explained in Hawkins 
et al (2003), an extensive analysis of redshift distortions in 
the 2dF catalog. 



see |http://lambda.gsfc. nasa.gov/product/map/dr3/parameters.cfm] 



To estimate the correlation ^(cr, vr), we use the ^ estima- 
tor of Landy and Szalay (1993) as in Ea. (|A5|l . with a random 
catalog Nr = 20 times denser than the SDSS catalog. The 
random catalog has the same redshift (radial) distribution 
as the data, but smoothed with a bin dz — 0.01 to avoid the 
elimination of intrinsic correlations in the data. Randoms 
also have the same mask. We count the pairs in bins of sepa- 
ration along the line-of-sight (LOS), tt, and across the sky, a. 
The LOS distance tt is just the difference between the radial 
comoving distances in the pair. The perpendicular distance 
between the two particles corresponds approximately to the 
mean redshift. It is exactly defined here as cr = y/s^ — tt^, 
where s is the distance between the particles. We use the 
wide-angle approximation, as if we had the catalog at an 
infinite distance, which is accurate until the angle that sep- 
arates the galaxy pair in the sky is larger than 15 degrees for 
the quadrupole and about 10 degress for the ^(a, tt) (Szapudi 
2004; Matsubara 2000). This later restriction corresponds 
to scales a > 80Mpc/h for our mean catalog. 

In Fig|4]we can see our £,{a, n) estimation for LRG DR6 
catalog and for MICE simulations (with a linear bias b=2 to 
see similarities visually). Here we can clearly see the FOG 
at small a as we have increased the pixel resolution when 
we approach the central part of the image (pixel size varies 
from 0.2 Mpc/h in the center to 10 Mpc/h at large scales). 
As expected, simulations have less noise because we have 
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Figure 4. 2-point anisotropic redshift correlation function ^{a, tt) 
for LRG galaxies in DR6 catalog (top panel) and mock MICE 
simulations (bottom panel, with a linear bias b=2, in order to 
be similar in amplitude to real LRG). Contours are -0.5 to -0.004 
with logarithmic bin of 0.4, 0. (over-plotted as a line), 0.003 to 40. 
with log bin=0.4. We see how the real LRG have a less flattened 
shape around the center, since /3 is smaller because of higher bias 



done the mean over 216 mocks. Moreover, the distortion pa- 
rameter /3 is higher in the simulations than in data because 
of the high bias of LRG, as can be seen in the shape of the 
2-point correlation function, which gets more flattened as /3 
increases. We see clearly the baryonic peak, as a ring in the 
correlation function ^{a, n). This feature is studied in detail 
in Paper IV of this series. In next sections we will analyze 
all the information hidden in this figure, using the analytical 
form of the error (see i)A2[) . 
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Figure 5. Q(s) (points with errors) and best fit asymptotic value 
of /3 at large scales (Eq. l(l8}) translated to the quadrupole (red 
dotted) 
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Figure 6. Top panel: Contour for (7„ and /3 obtained from the 
quadrupole Q(s) for the slice z=0. 15-0.47 at distances between 5- 
60Mpc/h. Bottom panel: Measured Q(s) (dots with error in gray) 
and overplotted the best fit (solid red) 
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4.1 Quadrupole and 13 estimation 

We use the normalized quadrupole Q(s) to obtain (5 — a„ 
(extensive explanation can be found in iiA5|) . But first, we 
can measure /3 using the asymptotic large scale value in 
the quadrupole where there is no dependence on cr„ (see 
Eq.llll)). We find /3 = 0.34±0.06 in the range 40-80Mpc//i; 
P = 0.32 ± 0.08 for 50 - 80Mpc/h and /3 = 0.34 ± 0.05 
for 40 — lOOMpc/h. In Fig. [5] we see the quadrupole with 
jackknife errors at large scales and the error obtained in /3, 
translated to the quadrupole as a band of red dots. 

The quadrupole only depends strongly on /3 and o"„. 
It does not depend on linear bias because it cancels out in 
the ratio, and it does not depend much on the shape of 
the 2-point correlation function (f2m and other parameters) 
or in the non-linear bias for small scales. We have tried to 
fit /3 and (j„ using all the scales in the quadrupole, fixing 
Qm to 0.25 (which in our model means fixing the shape of 
the real-space correlation function) and we use a power law 
form for the non-linear bias. When we change the shape of 
the ^(r) in the model (that is Qm for large scales and non- 
linear bias for small scales) , we obtain the same contours for 
P — Ov, SO we arrive to the conclusion that the normalized 
quadrupole is a good estimator to find /3 separately from 
the other parameters, which are degenerate with them in 
the ^(tt, a) plane. 

We can obtain better errors in /3 by doing a joint fit 
with (Jv and using smaller scales, as we do with simulations 
in Appendix A. We will fit the quadrupole above 5Mpc, 
because this seems a reliable scale according to simulations 
and results do not change much when using a cut on larger 
scales. Below this minimum scale, a„ starts increasing, and 
even when this has only a small effect in the quadrupole at 
larger scales, it can bias the value of 13 when doing a fit to 
all scales. 

In Fig|6]we show the fit as Ax'^ contours in the (3 — a^ 
plane (top panel) and the best fit to Q(s) for the mean slice 
(bottom panel). We have fitted /3 — CTi, for different slices 
in redshift. First, we divide the catalog in 3 redshift slices: 
z=0. 15-0.3, z=0.3-0.4, z=0.4-0.47. And then, we divide it in 
2 redshift slices: z=0.15-0.34, z=0. 34-0.47. The fitted values 
cr„ and 13 are similar in all the redshift slices (see table [TJ . 
This is a bit surprising because we expect / = Vlrniz)^'^^ to 
increeise slightly with z, for this redshift range. The similar- 
ity in the value of the distortion parameter (3 means that 
the bias must also increase as a function of redshift, roughly 
in the same way as / — Q.rn{zf'^^ so that the effect cancels 
out (see Eq.Q). This is not totally surprising, because we 
expect bias h to scale with redshift as 1/D{z), the inverse 
of the growth factor, D{z), which is proportional to /, in 
order to have stable clustering, ie b{z)D(z) constant along 
redshift (Peebles 1980). 



Sample 


/3 


cr^ (km/s) 


2=0.15-0.47 


0.310-0.375 


365-415 


z=0.15-0.34 
z=0.34-0.47 


0.280-0.365 
0.305-0.405 


320-410 
345-420 


z=0.15-0.30 
z=0.30-0.40 
z=0.40-0.47 


0.280-0.395 
0.285-0.365 
0.270-0.415 


305-435 
335-390 
305-420 



Table 1. Marginalized values for fS and (7„ to 1-a errors for each 
sample in redshift obtained from the study of the quadrupole 
Q(s). We can see that /3 is nearly constant along the different 
redshift slices 
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Figure 7. Real-space correlation function ^(r) (blue dots) and 
monopole in redshift space ^(s) for the main slice (orange dots). 
We see clearly how the redshift-space correlation function is the 
real-space correlation function biased by a constant factor that 
represents gravitational infall (dependent on /3 in Kaiser approx- 
imation) at scales above 4Mpc/h. However, for small scales, the 
redshift-space ^(s) is strongly suppressed compared to ^(r) due 
to random peculiar velocities. 



4.2 Fitting ^(r) 

We have calculated the real-space correlation function inte- 
grating through the LOS direction fEg. pTjl '). This method 
is complementary to the estimation from angular correlation 
(Baugh 1996). In Fig[7]we compare ^(r) (in blue) with the 
monopole ^(s) (orange) for the main slice in redshift. The 
difference at large scales (r > 5Mpc/h) is a constant value, 
at least until r < 30Mpc/h where we can trust the recovered 
(,{r). Note that, although we expect the recovered ^(r) to be 
systematically biased at large scales (eg see Fi^TJ, we can 
still see the baryonic peak at ^(r), in real space! So this is a 
strong feature in the LRG data. 

The large differences at small scales are due to random 
velocities, which we will study in more detail in Paper II of 
this series. 

In Fig|8]we show the ratio £,{s)/^{r) and we have over- 
plotted the constant value in the Kaiser approximation (Ea. (|23[l ). 
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Figure 8. We plot here the ratio 5(s)/5(r), and its prediction for 
large scales ('Eq. ||23|I . solid red line) using /3 = 0.34 obtained from 
the quadrupole Q(s) 



using the value of /3 obtained with Q(s), ie /3 = 0.34. This 
gives a good fit but note how it would be quite difficult to 
estimate /3 from this ratio given how noisy the data becomes 
when we reach an asymptotic constant value in C(s)/C(^)- 



4.3 Fitting large scales in ^{n,a) 

We obtain flm — Amp from the anisotropic redshift-space 
correlation function ^(cr, tt) using the full covariance matrix 
from the MICE group mocks (more details in i iA5|) . Here 
Amp refers to the factor b{z)as- The growth factor D{z), 
which is a function of Qrn, has been included in the model 
of ^(tt, a). We do not find any difference when fitting ^{a, n) 
with models with non-linear bias (as used in Paper-II) or 
linear bias. This is because we restrict our analysis to scales 
from 20Mpc/h to 60Mpc/h and for angles above 30-40 deg 
away from the line-of-sight in the a — tt plane. 

In Fig[9]we show constraints on Qrn. — Amp for the slice 
z=0.15-0.34. We have used the j3 — Ov obtained with the 
quadrupole Q(s) fit to marginalize the space of parameters 
Q.m — Amp — /? — o-„ in the ^{ti,(t) fit, although this only 
makes a slight difference. The values of /3 and (j„ are strongly 
degenerate in the range where we fit large scales, but cr„ is 
not strongly degenerate with the other parameters, since we 
are located away from the fingers of God, where we have the 
strongest effect of a„. When we marginalize, we assume that 
the likelihood scales as exp(— x^/2). 

We have done the same analysis for different redshift 
slices. We have summarized the results in table [2] where we 
have annotated the marginalized I-ct errors for the ampli- 
tude Amp=6a8 and for Q.m- Note that we use a fixed value 
of ft = 0.72 to model the shape the linear power spectrum. 

If we try to fit to even larger scales, we always obtain a 
slightly biased low Q.m compared to the one we obtain in the 
scales before the acoustic peak, probably because wide angle 







0.20 0.22 0.24 0.26 



0.28 



0.30 



Figure 9. Contours for Hm and Amp = berg, from a fit to g(7r, cr) 
on large scales, for the slice 2=0.15-0.34. We have marginalized 
for P — Oy, using priors from the quadrupole Q(s). Solid lines are 
1 — cr, 2 — cr and 3 — cr (1 dof), and dotted lines 1 — cr and 2 — cr 
(2 dof) 



Sample 


Amp = bag 


iim 


b 


f(iim) 


z=0.15-0.47 


1.47-1.65 


0.225-0.265 


1.73-1.94 


0.54-0.73 


z=0.15-0.34 
z=0.34-0.47 


1.45-1.62 
1.55-1.82 


0.230-0.275 
0.215-0.285 


1.71-1.91 
1.82-2.14 


0.48-0.70 
0.56-0.87 


z=0.15-0.30 
z=0.30-0.40 
z=0.40-0.47 


1.45-1.80 
1.42-1.60 
1.60-2.00 


0.240-0.320 
0.210-0.260 
0.195-0.305 


1.71-2.11 
1.67-1.88 
1.88-2.35 


0.48-0.83 
0.48-0.69 
0.51-0.98 



Table 2. Marginalized 1-sigma intervals for Amp = berg, Qm, b 
and f(Qrn) for each redshift sample. Here for b we used the best 
fit values of erg = 0.85 for the whole sample. 



effects that we are not taking into account in the modeliza- 
tion (see i]B5[) but possibly also because non-linear effects 
on the BAO peak and large sampling errors. We have also 
seen this effect in the simulations, where we have checked 
that the best region to obtain parameters is at intermediate 
to large scales, where we do our fit. 

Note how the best fit values of Amp = bag seem to 
change from sample to sample. This could be due to bias, 
which is both a function of z and luminosity. The values 
of Qrn agree within 1 — a for 2 degrees of freedom (dotted 
lines) . 

In Fig llOl we have plotted the £,{a, n) with the best 
model overplotted in solid lines. We can see how this simple 
Kaiser model can explain most features in the observations. 
There is a very good agreement on the region of negative 
correlation (in blue) which is a very good tracer of Qm and 
Qb (see Paper IV). All these features are very significant 
given the errors and their coherence over large regions. We 
can clearly see the FOG at small a as we have increased 
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Figure 10. Best fit (solid fines) to the 2-point anisotropic redshift 
correlation function ^(cr, tt) for LRG galaxies (colors). Colors are 
as in Fig |A12l l. The agreement with the model is excellent at all 
levels. 



the pixel resolution when we approach the central part of 
the image (pixel size varies from 0.2 Mpc/h in the center 
to 5 Mpc/h at large scales). The baryonic peak can also be 
seen in both the data and in the models. The significance of 
this is studied in detail in Paper IV, while the smaller scales 
and the prominent fingers of God are presented in Paper II 
(Cabre and Gaztanaga 2008). 



4.4 The value of ag 

We first try to obtain a fit to the parameter ag, which we can 
separate from the bias b{z) thanks to redshift distortions, 
following Eq. ^ and Eq. Q . We use our previous estimation 
of Clm — Amp from large scales (ie Table [U and the value 
of f3 from Q(s) (ie in Table [!}, using 7 = 0.55, for standard 
gravity. As Amp = b{z)as, we obtain as from: 

PAmp 

where Qrn{z) is given in Eq.([5ll. We also assume a fiat 
universe (WMAP results motivated) with a constant dark 
energy equation of state characterized by ui = — 1 (ie as in 
the cosmological constant model). 

In the top panel of Fig llll we show ag , marginalized over 
/3, Amp and il™ in each sample. We also show results using 
a fixed value of fim = 0.25 for all samples (bottom panel). 
There seems to be a lower ag detection for the middle slice 
of z=0.3-0.4, but it is consistent with the others at 2 — cr 
level. In the same figures we also show the estimated bias 



2.5 - 



1.0 



0.5 



m 


m 


WMM/m 


^^: 


'W' 


« 


mi//mi^/^^^ 


wMff//r 








~ 




- 



0.15 0.20 0.35 0.30 0.35 0.40 0.45 
redshift 



3.0 [ 



S 2.0- 



1.5 - 





0.5 L 



0.15 0.20 0.25 0.30 0.35 0.40 0.45 
redshift 



Figure 11. Estimation of linear bias, b(z), and ag for each slice 
in redshift. Gray: All the catalog, Red: z=0. 15-0.34, z=0.34-0.47. 
Blue: 2=0.15-0.3, z=0.3-0.4, z=0.4-0.47. In the top panel we have 
marginalized the values over the values of /3, Q.m and Amp = 
b{z)(Ts in each of the 3 samples shown. In the bottom panel we 
used a fixed best fit value of Qrn, = 0.25 and show results for all 
subsamples. 



b{z) from the amplitude Am,p and the values erg. For a fixed 
Qrn, the linear bias b{z) seems higher as we move to higher 
redshifts, and it is consistent with previous results found also 
with LRG (Zehavi et al 2005; Tegmark et al 2006; Okumura 
et al 2008). Note that the luminosity of the galaxies also 
increases with redshift and this is probably the main reason 
for the increase in b (more luminous galaxies trace higher 
density peaks and therefore are more biased). There is one 
exception to this tendency in slice z = 0.30 — 0.40. The 
values of b for the best fit ag = 0.85 are shown in Table [2l 
were we can see how 6 is lower for z = 0.30 — 0.40. This is 
due to the LRG selection explained in |3l which includes a 
population of less luminous galaxies Mr > —21.5 at z ~ 0.35 
(see Fig |B12|l resulting in a lower bias. 

Once we know the bias we can use the values of /3 to 
estimate the linear velocity growth function f{Q,rn) for each 
slice, which are also shown in Table [2] 
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4.5 Modified gravity 

We can also use redshift distortions to constrain modified 
gravity. For standard gravity, the linear theory growth fac- 
tor, D{a) — 5/5(0), depends purely on the expansion history 
H[a), w{a), ilm{a.)- Any discrepancy found between the ob- 
served growth factor and predictions based on the expansion 
history can be used to test the gravity. 

The idea of using the growth of structure to test grav- 
ity has a long history and dates back to Brans-Dicke (BD) 
model (see Brans 2005 for a historical review) . More recently, 
after cosmic acceleration, the BD and variations have been 
revisited to explore structure formation outside general rela- 
tivity (Gaztanaga and Lobo 2001; Lue et al 2004; Lue et al 
2004; Acquaviva et al 2008). Zhang et al (2007) propose 
a test to discriminate between different models for grav- 
ity on cosmological scales based on lensing. Nesseris and 
Perivolaropoulos (2008) have compiled a data set of various 
observations of the growth rate / at different redshifts, ob- 
tained through redshift distortions or indirectly through the 
rms mass fluctuation as{z) inferred from Ly — a. These data 
points are used to constrain and parametrize the linear per- 
turbation growth rate /. Wang (2007) does a prediction of 
the characteristics that a survey must accomplish to be able 
to rule out the DGP gravity model (an extra-dimensional 
modification of gravity), where the idea is to calculate H(z) 
from the baryon acoustic peak and f(z) from redshift distor- 
tions. Guzzo et al (2008) test the nature of cosmic accelera- 
tion using galaxy redshift distortions at z=0.8, obtaining /, 
but errors are still too high to distinguish between different 
theories. Acquaviva et al (2008) have recently done a new 
compilation of results. See Bertschinger and Zukin (2008) 
for a theoretical approach to modified gravity. 

If we fix a value for erg , which can also be known from 
other observations, we can assume that the changes in the 
^(tt, a) amplitude with redshift could be explained by changes 
in the growth factor D{z) due to a different law of gravity at 
cosmological scales (Linder 2005). This can be represented 
by the growth index 7. Both f{z) and the growth factor 
D{z) changes with 7. We have plotted in Fig |12l how / and 
D change with 7 for z=0.34 (mean redshift in LRG) and 
^m = 0.25. The measured ^(vr, a) amplitude on large scales 
depends linearly on the product asb{z)D{z). If we now fix as 
we can have an estimation of b{z)D{z) as well as our sepa- 
rate estimation of (3, from the quadrupole. We therefore can 
use this combined data to produce an estimate of f{z)D{z). 

We have assumed that the observations can be explained 
by varying 7, for a fixed dark energy equation of state w. 
This is a good approximation because w only depends slightly 
on 7 (Linder 2007). We show in Fig |13l our estimation of 
1 — (T errors for the growth index 7 once we fix a8=0.7, 0.8, 
0.9 and marginalized over /3, Amp and fim- As shown in 
FigdU the product fi'y)D{j), at a given redshift, decreases 
with 7. If we change the value ag to higher values, the fac- 
tor P b{zsiice) D{zaiice) obtained from observations will be 
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Figure 12. We see here the change of / and D with 7 and the 
factors involved in redshift distortions and ISW effect. We have 
fixed Qrn = 0.25 and z = 0.34, as in the main LRG catalog 



lower, and this needs to be compensated by reducing / D, 
thus increasing 7. In next section ( ^4.6p we will show how 
the argument goes in the opposite way when we work with 
ISW effect. At 2 — (T and for all the redshift slices, 7 is con- 
sistent with a standard gravity, except for ag = 0.7 where 
we need 3 — cr for the last slice, favoring a ag clearly higher 
than 0.7, which is in agreement with recent observations. 

Following this argument, we can calculate which is the 
range of ag that is consistent with standard gravity, 7 = 
0.55. We find 0.80 < erg < 0.92 at 1-cr. In the modified grav- 
ity DGP model, erg needs to be higher since 7 is also higher 
(7 — 0.68). We could say that DGP model is inconsistent 
with data if ag < 0.84. These results are in agreement with 
the recent study of Yamamoto et al (2008), which appeared 
as we were writing this paper. 



4.6 Angular correlation and cross-correlation with 
WMAP 

We can also learn about the growth of fiuctuations using the 
ISW effect. This method is independent to the one explained 
before with redshift distortions (see Cabre et al 2007 and 
references therein for an explanation of the ISW effect in 
the context of SDSS galaxies). In this section we explore the 
angular correlation function and also the ISW effect through 
the cross-correlation between galaxies and fiuctuations of 
temperature in WMAP. Fig |14l shows the results for wgg 
and wtg at the redshift slice z=0.15-0.34. This is the same 
data that we have used before, but we now look at angular 
clustering. The angular auto-correlation function GG scales 
as: 



WGG oc as(l)G{z) b{z) D{z) , 



(25) 

while the cross-correlation function between galaxies 
and CMB temperature fiuctuations wtg is proportional to 
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Figure 13. Growth index 7 for different redshift slices as in top 
panel of Fig llll when we fix a^ = 0.7,0.8,0.9 (from top to bot- 
tom). As we increase trg we find that 7 also wants to be higher, 
in favor of modified gravity, such as DGP, where 7 = 0.68 (in 
standard gravity 7 = 0.55) 



WTG OC oUG{z)h{z)D{z)'^^^^'^'''^ 



dz 



where 



d[D{z)/a] 
dz 



^D{z){l-f) 



(26) 



(27) 



and <j>g{z) is the galaxy selection function. Both wgg and 
wtg are proportional to al because this factor conies from 
the normalization of the power spectrum, but wgg is pro- 
portional to {(j)G{z)h{z)D{z)) while wtg is proportional 
to {(j)G{z)h{z)D{z)) (from the clustering of galaxies) and 



— — ]^^ (from the evolution of gravitational potentials). 
Thus, this allows for a separate estimation of these quan- 
tities (Fosalba et al 2003). 

We find that the measured signal of wtg in Fig |14l is 
higher than expected, a clear tendency that has been seen 
before (see Gaztanaga etal 2006, Giannantonio et al 2008 for 
a compilation of ISW observations). The high signal wtg 
could be due to: sampling variance, higher ug, lower n^, 
non-linear effects, bieis between matter and galaxies differ- 
ent from the one obtained from galaxies-galaxies, non-linear 
bias, different form of dark energy as w > —1 (see Cabre 
et al 2006 for some hints in this direction) , modified gravity 
at cosmological scales, or non-linear magnification (linear 
magnification is not expected to affect ISW at low redshifts, 
Loverde et al 2007). 

The signal to noise is not very high so we can not obtain 
tight constraints, but we can explore what could be creating 
this high signal. We study two reasons here: a change in as or 
a change in the growth index 7. We have also studied these 
two parameters in the section above of redshift distortions 
and we want to see if results are compatible. 

We can break the degeneracy between b{z) and ag in 
the auto-correlation function wgg, which is proportional to 
6^(t|, by combining the result with wtg, which goes as fccrg. 
We will assume that b{z)D{z) = b{Zsiice)D{zsUce) is con- 
stant through all the redshift slice, to be consistent with the 
previous section of redshift distortions. 

We fix the shape {Qm = 0.25 and fiat universe) of wgg 
and use the amplitude to find the factor b{zsUce)D{zstice)o-s- 
This should be equal to the amplitude that we found in 
previous sections, when analyzing 3D redshift distortions, 
since we are working with the same LRG galaxy samples. 
We find that this is the case: both measurements of the 
amplitude are consistent within the errors. There is an extra 
power at large scales (> 4 degrees) as in the 3D case that 
can be due to sampling variance. This is explained with more 
detail in ^B] We next fix the dark energy equation of state 
parameter to w — —1 and 7 = 0.55 to standard gravity. 
We now try to explain the observed amplitude of wtg by 
just changing as- We therefore break the degeneracy of erg 
with bias b{z) in a way that is completely independent from 
the one used in the previous section, based on 3D redshift 
distortions. In bottom panel of Fig |14l we have compared 
Wtg to the best model (in red) which corresponds to ag = 
1.2. This value is quite high but the errors in wtg are also 
large, so that at 1-a (or 2 — a for some of the redshift slices), 
the values of as are consistent with our previous estimation. 

We next see how we could explain the observed high 
ISW signal if it is due to a modification of gravity. We fix 
ug = 0.8 and vary 7. We assume that w{z) changes only 
slightly with 7 and that almost all the variation with 7 in 
ISW comes from the factor D(l - /). Both D and (1 - /) 
grow with 7, so that D{1 — f) also grows with 7 (see Fig |12p . 
Thus if we fix as to a low value we need larger values of 7 
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Figure 14. Top panel: wqg of LRGs (solid with errors) and best 
linear model (dotted line) for the slice in redshift z=0.15-0.34. 
Bottom panel: Wtg (solid black with errors) and best model when 
fitting (Tg, for a (T8 = 1.2 and standard gravity, 7 = 0.55 (solid 
red line). 



from the ISW effect, but lower values of 7 from redshift 
space distortions. Thus, this seems to be a promising test 
because both ways to obtain the growth histories seem to 
be constrained by data in opposite directions (ie see Fig|12p. 
and we therefore can use these test to break degeneracies in 
future surveys with better ISW signal-to-noise. For the slice 
2=0.15-0.34 and for & as, = 0.8, we find that 7 needs to be 
7 = 0.8, but errors are in this case really big, fully consistent 
with standard gravity. 

From these observations, we conclude a preference for 
either a higher cts or equivalently, a higher 7 than in stan- 
dard gravity. But this is only a 2-sigma effect. High power 
could also be due to cosmic variance. Fundamentally, we 
want to remark that ISW can provide independent and com- 
plementary information to the growth parameter 7. This is- 
sue will be resolved with future surveys. 



4.7 Different redshift slices 

We next look at the differences between the redshift slices 
in the monopole (Fig llSp . in the projected correlation func- 
tion H((t) (Fig|16[l. in the real-space correlation function ^(r) 
(FigHlll and in the ratio C(s)/^(r) (Fig fTS]) . 

As we have seen in section ^4.11 13 is similar for all the 
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Figure 15. Comparison between ^(s) in redshift space for differ- 
ent slices in redshift: All (black), z=0. 15-0.3 (solid blue), z=0.3- 
0.4 (dashed blue), z=0.4-0.47 (dashed-dotted blue); z=0.15-0.34 
(solid red), z=0.34-0.47 (dashed red) 
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Figure 16. Comparison between the projected correlation func- 
tion H(cr) for different slices in redshift. All: black, z=0.15-0.3 
(solid blue), z=0.3-0.4 (dashed blue), z=0.4-0.47 (dashed-dotted 
blue); z=0.15-0.34 (solid red), z=0.34-0.47 (dashed red) 



redshift slices. We can also see this in Fig llSI which shows 
that the ratio i,[s)/(,{r) is quite similar at large scales, indi- 
cating similar values of 13 (Ea. (|23[) ') in all samples. In some of 
the cases, there seems to be a turning down of the ^(s)/^(r) 
ratio at scales larger than 30 Mpc/h. This is due to the bias 
in the recovered ^(r) which does not work well for scales 
larger than 30Mpc/h (ie see FiglTJ. 

The monopole is approximately a measure of the shape 
of real-space correlation function for large scales, but scaled 
up by a function of 13 similar for all the slices. Looking at 
the monopole (Fig llSf) and also at the projected correlation 
function (Fig llGp and at the real-space correlation function 
(Fig llTp we see that all the slices except from the further one 
(blue dash-dot) lay in the same line, meaning that D{z)b{z) 
is almost constant with redshift, what is sometimes called 
stable clustering. 

In Fig |17l we also see that the change of the slope at 
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Figure 17. Comparison between ^(r) in real space for differ- 
ent slices in redshift. All: black, z=0.15-0.3 (solid blue), z=0.3- 
0.4 (dashed blue), z=0.4-0.47 (dashed-dotted blue); z=0. 15-0.34 
(solid red), z=0.34-0.47 (dashed red) 
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Figure 18. Comparison between i;{s)/^{r) for different slices in 
redshift. All: black, z=0.15-0.3 (solid blue), z=0.3-0.4 (dashed 
blue), z=0.4-0.47 (dashed-dotted blue); z=0.15-0.34 (solid red), 
z=0.34-0.47 (dashed red). The constant solid black line shows the 
Kaiser expression for ^{r)/^{s) at large scales (as a function of /3, 
see Eq.llllll) for /3 = 0.34 



small scales moves to larger scales as we more to more dis- 
tant slices. We think that this effect is probably not physical 
but due to geometry since it corresponds to a fixed angular 
scale, ie the mask or the fiber collision problem. 

Fig |19l shows a zoom over the BAO peak region, as mea- 
sured in the monopole. This peak is roughly in agreement 
with predictions and is dominated by modes in the perpen- 
dicular direction. The BAO peak is detected in all slices 
and our best fit model (shown as a continuous line in the 
Figure) seems a bit lower than the data. This can be fixed 
by exploring the full parameter space of cosmological mod- 
els and will be done in Paper IV of this series. This is an 
important test because it reproduces previous results and 
indicates that systematic effects can not be too important 
in our analysis on scales of 100 Mpc/h or smaller. On larger 
scales, some of the slices show more power than predicted 



Figure 19. Comparison between the BAO peak in ^{s) for dif- 
ferent slices in redshift. All: black dots with errors, z=0.15-0.3 
(solid blue), z=0.3-0.4 (dashed blue), z=0.4-0.47 (dashed-dotted 
blue); z=0.15-0.34 (solid red), z=0.34-0.47 (dashed red). These 
measurements are compared to the best fit linear model for the 
principal redshift slice (solid line) 



by models. This is compatible with the errors, but we won- 
der if there could be some hidden systematics here. We will 
explore this further in the Appendix B which is dedicated 
to the study of systematic effects. 



5 DISCUSSION AND CONCLUSION 

In this paper we work with luminous red galaxies (LRGs) 
from the catalog SDSS DR6 to obtain cosmological parame- 
ters from large scales. These are good galaxies to study this 
range of scales since they trace a big volume of the universe. 
We use a simple approximation to model the redshift space 
distortions (see ^, which work incredibly well with LRGs 
on large scales, as explained in this paper, but also on small 
scales, shown in Paper II (Cabre and Gaztanaga 2008) of 
this series. 

We have used one of the largest N-body simulation run 
to date (with box size of L=7680Mpc/h with 2048^ parti- 
cles) to create realistic mock LRG galaxy catalogs to analyze 
the data (Appendix A). We have checked our results to some 
extent with smaller simulations. This is a crucial step given 
the complications in the data and the accuracy needed for its 
interpretation. We have validated all the methods that we 
use during the paper with credible mock simulations. We 
also test the validity of the models and the errors: Monte 
Carlo, jackknife and a theoretical approach for ^(a, vr). We 
find that JK errors does not work properly for large perpen- 
dicular distance a. 

There are usually two important sources of statistical 
errors in clustering analysis: Shot-noise error, due to the fi- 
nite number of galaxy pairs, and sampling errors, due to 
the finite amount of volume in the survey. We have shown 
with our mock simulations that the LRG samples in SDSS 
DR6 are dominated by the shot-noise term in the 2-point 
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5(a, tt). Usually, this term is assumed to follow the Poisson 
distribution, with an error given by the inverse of the square 
root of the number of measurements (ie see Ea. (|A6|) ). This 
is in fact the case to a good accuracy in our dark matter 
mocks. But this is not the case for groups of dark matter 
particles (halos) which are selected in a way that, by con- 
struction, the detection of a group creates a large exclud- 
ing region around the group position. This produces what 
we have called super-Poisson errors, which are about 40% 
larger than Poisson errors. Jack-knife errors are shown to be 
good at detecting this additional shot-noise component and 
we find that this effect is also present in the real LRG data. 
This is interpreted here as direct evidence that the domi- 
nant fraction of the LRG galaxies populates separate dark 
matter halos. But this does not exclude a minor fraction of 
them to be found with other LRGs within the same halo. In 
fact, in Paper III, (Gaztanaga et al 2008), of this series we 
find evidence for a large quadratic bias which indicates that 
this is the case (Scoccimarro et al 2001). 

We conclude our discussion on errors by emphasizing 
the importance of using realistic mocks to assess our results. 
We have demonstrated here that dark matter mocks are not 
good enough. Diagonal, sampling errors can in principle be 
scaled from dark matter to real data because the relative er- 
ror only depends in the sampling volume. But this is not the 
case for the shot-noise or for the covariance. The difference is 
important. Errors and covariances are underestimated with 
dark matter mocks and results will be wrongly interpreted. 

We have calculated the distortion parameter /? = [0.310— 
0.375] using the normalized quadrupole Q(s), which also pro- 
vides us of an effective value for the dispersion of random 
peculiar velocities (t„ = [365 — 415] km/s (see Table [l] for 
more slices in redshift) . The great advantage of the normal- 
ized quadrupole introduced by Hamilton (1992) is that it 
can measure the squashing produced by bulk galaxy mo- 
tions with independence of the overall amplitude or shape 
of the galaxy 2-point correlations. We have also checked if 
this value for /3 gives the correct ratio £,{s)/^{r) on large 
scales. We have shown that with actual data, we can only 
recover the real-space correlation function ^(r) for scales 
smaller than 30Mpc/h, and we see that the ratio between 
the redshift-space monopole ^(s) and real-space correlation 
function ^(r), ^{s)/£,{r), can be described by the approxi- 
mated Kaiser value at large scales, with the same value of j3 
found in the quadrupole. 

We obtain a value for flm from the shape of the redshift- 
space correlation function in the anisotropic plane tt — u at 
large scales. The value is Q,n = [0.225 — 0.265] for the mean 
slice (see Table [2] for other redshift slices). We also obtain 
a value for the amplitude bag in the linear regime. We can 
therefore break the degeneracy between as and b thanks to 
redshift distortions through the parameter /3, which results 
in as = [0.79 - 0.91] and b around 2. 

We next look at the growth history and possible mod- 



ifications of gravity through the parameter 7 which is 0.55 
for standard gravity and 0.68 for the modified gravity DGP 
model. We find consistence with standard gravity for 0.8 < 
as < 0.92 at 1-a. DGP model is inconsistent with our results 
ff as < 0.84. 

We have also cross-correlated LRGs with WMAP in or- 
der to investigate the ISW effect, obtaining a high signal, as 
in other studies, compared to current A.CDM model. The 
degeneracy 6 — erg and the growth history can also be studied 
using the ISW effect which is independent from the study of 
redshift distortions. The cross-correlation indicates a prefer- 
ence for higher as- 

In the Appendix B we look for systematic errors, and 
arrive to the conclusion that our results are robust in the 
scales we are doing the analysis. 
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APPENDIX A: SIMULATIONS AND ERRORS 

We first present the simulations used in this paper. Next we 
show how simulations have been used to provide an error 
estimate for the different statistics used in this paper. Fi- 
nally we use the simulations to validate the methods that 
we will apply to real data. We will conclude that the meth- 
ods and errors that we use to obtain parameters from LRG 
data are all well validated by simulations. We follow similar 
methodology and notation as in Cabre et al (2006). 



Al Simulations 

We have used different comoving outputs (between z=0.0 
and z=0.5) of a MICE simulation, run in the supercom- 
puter Mare Nostrum in Barcelona by the MICE consortium 
(www.ice.cat/mice) in order to study the errors and vali- 
date the models used (Fosalba et al 2007). The simula- 
tion contains 2048^ dark matter particles, in a cube of side 
7680Mpc/h, ^M = 0.25, fib = 0.044, ag = 0.8, n, = 0.95 
and h — 0.7, with a flat cosmology. We have divided this 
big cube in 3^ cubes of side 2 x 1275Mpc/h, and taking 
the center of these secondary cubes as the observation point 
(as if we were at z=0), we apply the selection function of 
LRG, which extends to z=0.47 (r=1275Mpc/h). We can ob- 
tain 8 octants from the secondary sphere included in the 
cube, so at the end we have 8 mock LRG catalogs from each 
secondary cube, which have the same density per pixel as 
LRG in order to have the same level of shot noise, and the 
area is slightly smaller (LRG occupies 1/7 of the sky with a 
different shape). The final number of M independent mock 
catalogs is 216 (27x8). 

We use both dark matter particles and groups in this 
simulation. By definition, there is no bias in the dark mat- 

ter particles, so /3 = "^\^' = 0.62, where Qm(z) — 

nm(l + zf/{nm{l + z)'-' + 1- Qm) and b = 1. To simulate 
biasing we select groups of particles using friend-of-friends 
with linking scale of 0.20 (we have also used 0.16 with consis- 
tent results). At z = we find a total of 107 million groups 
with more than 5 particles (M > 1.87 x 10^^). When the 
number of particles in the group is larger than about 40 or 
so, they usually correspond to DM halos and have similar 
mass function (Evrard et al 2008) and bias. We find that 
50 — 80% of groups with 5-40 particles (respectively) cor- 
respond to DM halos. The rest of the groups are just high 
density regions, which in any case are also candidate sites 
for bias galaxy formation. We therefore use groups of differ- 
ent masses as biased tracers and use the mass threshold to 
choose the group subsample that better matches the number 
density and clustering amplitude in the real LRG data. We 
have also used a MICE simulations with 2048^ dark mat- 
ter particles, in a cube of side 3072Mpc/h (which we call 
MICE3072, same parameters as MICE7680) which has 15 
times better mass resolution to check for mass resolution ef- 
fects. We find very similar results in both cases to the extent 
that we can compare the smaller number statistics. When 
we select groups with M > 2 — 4 x 10^^ , both the clustering 
amplitude (6i ~ 1.9 — 2.2 for as — 0.8) and the number den- 
sity (n ~ 4— 6x 10^"") are similar to the real LRG galaxies in 
our SDSS sample (the range refiects the fact that the actual 
number depends on LRG sample used, ie redshift and selec- 
tion). We have explored a number of cases to make sure we 
understand how errors depend on number density and bias. 
We will focus in presenting results for 3 different cases: a) 
dark matter particles aX z = 0.3 {b = 1, P — 0.62), diluted 
to fit LRG number density b) groups with M > 2.2 x lO" 
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(6 = 1.9, /3 = 0.25) at z = 0, c) groups with M > 3.74 x 10^^ 
(b — 3, (3 = 0.23) at 2 = 0.5, which cover the range b, /3, z 
and densities in the real LRG in our samples. 

We generate the mock catalogs applying redshift distor- 
tions in the line-of-sight direction. 



-Vr/H{z)/a(z) 



(Al) 



We obtain the Monte Carlo (MC) error from the dis- 
persion of M independent realizations of our universe. For 
M reahzations, the MC covariance is 



C-i 



1 ^^ 



-m){m-m) 



(A2) 



where £,{i)^ is the measure in the k-th simulation (k=l,...M) 
and ^(i) is the mean over M realizations. The case i=j gives 
the diagonal error (variance). Typically, we need M = 100 
independent simulations for the diagonal error (in order to 
obtain 5% accuracy), and more for the covariance matrix, 
depending on the case. The MC error is an estimation of 
the true error, but it takes lots of computational time and 
it also requires simulations with the same particularities of 
the data analyzed. When we refer to off-diagonal elements 



of Cij we use the normalized covariance: 



Cij 



Ci- 



yCiiCjj 



(A3) 



which ranges from -1 to -fl. Values of \Cij\ < 0.2 usually 
have little impact on the error analysis. 

In the case of the jackknife (JK) error, we obtain the 
different realizations that we need to compute the error from 
the same data (or from a single realization in the case of sim- 
ulations). This is done by dividing the sample in M zones. 
Each JK subsample consists in the whole catalog except 
from one of these M zones. In this case, as the realizations 
are clearly not independent, we apply a multiplicative fac- 
tor of (A/ — 1) to the previous covariance to account for this 
effect. Fig lAll shows the octant catalog used for the galaxy 
mocks divided in 63 jackknifes zones, with similar area and 
shape; and the SDSS DR6 real catalog divided in 73 jack- 
knife zones. Note how the real mask is different than the 
mask in the mocks, but note that the JK regions are similar 
in shape. We have tested that this difference in mask does 
not make a large difference in the conclusions of our error 
analysis, so we find adequate to work with an octant which 
allows us to obtain a higher number of mocks, necessary for 
a good determination of the error. We have studied with a 
reduced set of mocks with the proper LRG mask that the 
general shape of the mask is not important for the large 
scale analysis as far as it remains compact. Small holes in 
the mask do not afi'ect the errors at large scales. In this 
way we can apply the conclusions obtained from our mocks 
to the real LRG data. We use the dark matter and group 



JK zones in DR6 catalogue 




JK zones in the mock catalogue 




Figure Al. Top panel: JK zones for the LRG catalog (equatorial 
coordinates). Bottom panel: the mock catalog with its JK zones 



mock simulations to probe the limit in which the JK errors 
represent a good approximation. We then use the JK error 
from the real LRG data when adequate. We study the er- 
rors in the redshift space correlation function ^(7r,a), the 
monopole ^(s) and the quadrupole Q{s). Finally, we use the 
simulations to validate the models that we will use with real 
data. 



A2 Errors in ^(a, vr) 

We study here the errors in the anisotropic redshift-space 
correlation function ^{a,n). The JK and MC errors only 
seem to agree at small perpendicular scales a < 20Mpc/h, 
which means that we can not use the JK error for large 
a values. This is illustrated in Fig |A2l The reason for this 
behaviour could partially be related to the fact that we have 
defined the JK zones in angular space, ie in the a direction. 
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Figure A2. We compare diflferent error estimators for tiie 2-point correlation function in redshift space ^(cr, vr) for MICE dark matter 
mocks (left panel) and group mocks with b = 1.9 (right panel). We fix the perpendicular distance a and move along the line-of-sight vr 
as indicated in the figure. The errors are: MC (solid line), theory (dashed line) with and without the signal part of the error, and JK 
(dotted line with errorbars). As we move to higher tr (lower amplitudes), JK starts to fail, but the analytical and MC errors always agree 
well. On the radial direction, top lines for a = 2.5Mpc/h, the JK and MC agree well, but theory model only works in the dark matter 
case. 



As we increase a the difTerent JK zones are less independent. 
This could explain why JK errors do not work so well. In the 
radial direction, the different JK zones are more independent 
from each other and the JK errors give a better agreement 
to MC errors. 



A2.1 Error model 

We have found a phenomenological model (we sometimes 
call it "theory" or "analytical" error) for the errors that 
match well the MC errors at large a, even at large scales. 
The advantage of using a model is that is smoother than 
MC errors and can also be used to calibrate JK errors. We 
will first test this model in different situations. 

Our "theory error" has a dominant shot-noise contribu- 
tion, which scales as 1/ ^number of pairs, and a part pro- 
portional to the signal. This works well in this case because 
our LRG sample is shot-noise dominated, specially at large 
scales, and we are able to separate both contributions, as 
if we were in a diagonal space. In general, the analytical 
derivation of the error is more complicated that this simple 
modeling. We propose the error to have the following form 
A^ — A£^shot~noise+A£^signai, with two arbitrary coefficients 
cinoise and Osignai , SO that: 



AC = a„ 



Mf 



I Oisignal S 



(A4) 



adding these two terms linearly fits the simulations better 
than in quadrature. For Poisson shot-noise we have by defi- 
nition that anoise = 1. This works very well for dark matter 
mocks, but we will show that groups do not follow the Pois- 
son distribution and Onoise = 1.4. We can associate Osignai 
to the number of independent nodes in the catalog and we 



therefore expect Usigmai to scale with the square root of the 
volume of the sample used. We also expect a^ignai to depend 
on the binning used. We use the ^ estimator of Landy and 
Szalay (1993) to calculate the correlation function, 



a<^,^) 



DP - 2DR + RR 
RR 



(A5) 



In the case of dark matter mocks, the error that comes 
from having a limited number of data-data (DD) pairs (Pois- 
son shot-noise) is: err{DD) = l/^RR/Nf^, and for DR 
pairs: err{DR) = 2/ ^/RRJNr, where the random catalog is 
Nr times denser than the data catalog. The error in random- 
random is insignificant, because we are using a denser ran- 
dom catalog. We add the different errors in quadrature to 
find: 



Ae 



Poisson — 



iV|. 



-4iVH 



RR 



(A6) 



We can reduce the error in ^(tt, a) by increasing the 
number of particles in the random catalog, modifying the 
error coming from DR, which is inversely proportional to 
the square root of the pairs data-random. The first part of 
the error is always the same, because it depends only on 
the data. For A^'^ = 10, as in the simulations, A^poisson = 
^UO./RR, so the DR part is 40% of the DD part. For Nr = 
20, as in the analysis of LRGs, A^poisson = \/4S07/RR, so 
the DR part is now 20% of the DD part. For a binning of 
5Mpc/h, we find asignai = 1/95 and for binning of IMpc/h, 
O-signai = 1/25. For (T > 20 Mpc/h, the signal part in ^(tt, a) 
is not very significant. This can be seen in Fig |A2l where the 
linear dashed line correspond to the shot-noise term and the 
curved dashed line (only important on scales tt < 20Mpc/h) 
also includes the signal part in Ea. (|A4p . For small a scales. 
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Figure A3. Diagonal error in ^(cr, n) in redsiiift space for MICE darli matter simulations (top panel) and groups with b = 1.9 (bottom 
panel) with contours § = 0.0035 — 0.05 (log increment=0.1) 



JK and MC coincide, so we could just use JK error from 
the actual LRG data. In practice we prefer to use our model 
error to the JK error because this avoids introducing noise 
in the error analysis, but we have checked that none of our 
results depends on this choice. 



A2.2 Super-Poisson errors 

In the case of groups, the Poisson model for the shot-noise 
does not work well. The reason for this is that groups are se- 
lected using friend-of-friends with linking length of r = 0.20 
times the interparticle separation. This means that groups 
create excluding regions where we can not randomly locate 
another group and therefore do not follow a Poisson distri- 
bution (Smith et al 2007). Shot-noise is larger in this case, 
which we call "super- Poisson". After exclusion, the fraction 
of available volume is (1 — r)'' ~ 0.5 so we expect the Poisson 
term to be corrected by the square root of this fraction, ie 
roughly a„oise — 1/(1 — r)'^'^ ~ 1.4. This is exactly what 
we find for our model in Eq. (|A4|l for group mocks. We have 
studied super-Poisson errors in two group simulations (z=0 
and z=0.5) with the same linking length r = 0.20, and ob- 
tain the same value of Onoise — 1.4, as expected, despite the 
large difference in the number density of groups. 



What about real LRG galaxies? Do they follow the Pois- 
son shot-noise or the groups shot-noise model? Real LRG 
data seems to follow more closely the group shot-noise as 
we find that we need cvnoise = 1.4 rather than anoise = 1 for 
Poisson term when we fit the JK errors in the real data at 
small a. As we have seen with mocks, the JK errors repro- 
duce very well the true errors for small a both the Poisson 
(dark matter) and super-Poisson (groups). The amplitude 
of the JK errors on small a is quite closer to group mocks 
than to dark matter mocks (this can be seen by comparing 
Fig lASI to Fig |A4l below). This by itself is a very interest- 
ing result because it clearly shows that LRG galaxies have a 
tendency for exclusion. This suggests that a dominant frac- 
tion of LRG populate separate dark matter halos (this has 
already been studied with the same conclusions using semi- 
analytical models, Almeida et al 2008, or with real LRGs, 
Wake et al 2008 and Zheng et al 2008). It also shows that us- 
ing dark matter mocks or a standard (Poisson) point process 
to simulate LRG will result in an important underestimation 
of errors. 
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Figure A4. Same as Fig |A3l for real LRG data. 



A 2. 3 Errors in LRG 



Fig |A3l shows the full 2D version of Fig |A2l We show here 
the differences in the diagonal error between MC, JK and 
the theory error in Ea. (IA4|) . when binning the correlation 
function with 5Mpc/h. MC error is calculated using MICE 
simulation and JK is the mean over all the JK that we have 
calculated in each mock. JK errors work well for small a < 
20Mpc/h but that become higher than MC when going to 
large a for the reasons explained above. Our "theory" error 
model form agrees with MC error at all scales, except in 
the radial direction, a = 2.5Mpc/h, for the case of group 
mocks (where we used a fixed Qnoiae = 1.4 at all a), which is 
slightly higher than the MC or the JK errors. The model for 
dark matter mocks, with anoise = 1 does well at all a. This 
can be seen in Fig |A2l Thus in the radial direction we need 
to change anoiae for groups to anoise = 1-2, a lower value 
than anoise ~ 1-4. The reason for this is most probably the 
random component of peculiar velocities in redshift space, 
which ramdomizes the radial positions of group and reduces 
the super-Poisson shot-noise contribution, which comes from 
self-exclusion in group positions. This is no problem when 
analyzing real data as we can fit a„oise to match the JK 
error at small values of a (with a different value at cr < 5 
than at 5 < CT < 20 Mpc/h) and then use our error model 
for large values of a where JK does not work well. We could 
also directly use the MC errors in the group mocks, which 
seem to agree quite well with real LRG data (ie JK errors 
are very similar in both cases, see Fig |A4p . 

In Fig |A4l we show the JK error obtained from the real 
LRG data using a random catalog 10 times denser than the 
data (as done with the mocks), and the analytical form with 
a fixed anoise = 1.4 and cxsignai ~ 1/95. As in the mocks, 
for large a the JK error is bigger than the analytical error, 
which should be more representative of the true (MC) error. 
At small scales the model follows the JK prediction, except 
for (7 < 5 Mpc/h where we need to change to Ofnoise = 
1.2 (the figure shows the model for anoisf^ = 1-4 at all a). 
Note the similarity between the errors in the data and in the 




30 30 

TT(Mpc/h) 



Figure A5. Normalized covariance of §(71, a) for groups in a bin 
centered at position (a = 32Mpc/h,-K = 28Mpc/h). Contours 
show fix values of the covariance away from the central position, 
with contours 0.5, 0.3, 0.2 and 0.1. 



group mocks in Fig lA3l For our final analysis of the data we 
use a random catalog 20 times denser to reduce the error 
contribution due to the shot-noise in the randoms. 

We have done a similar analysis with the L-BASICC 
halo simulations from Durham (Angulo et al 2008), which 
have approximately the same bias as real LRG in SDSS, 
and we find very good agreement with our fitting formula. 
If we use a random catalog 20 times denser, as we do for the 
analysis of the data, the error is lower, just as expected. 



A2.4 Covariance (,{tt, a) 

So far we have shown diagonal errors. The normalized covari- 
ance is in fact quite small for dark matter mocks. It has the 
same approximate shape and amplitude for all the points. 
The amplitude of the covariance is just given by the distance 
between bins: at a separation of 5Mpc/h the normalized co- 
variance is < 0.1. This is also true for group mocks in the 
radial direction. But the covariance is larger in the case of 
groups for a > 5Mpc/h. At a radial distance of 5Mpc/h 
between bins in ^(tt, a) the normalized covariance is about 
0.5 and dies below 0.2 for bins separated by more than 15 
Mpc/h. This is illustrated in Fig |A5l which just shows the 
covariance in one of the bins. Other cases are very similar. 
The inner contour, which corresponds to 0.5 for groups, is 
very similar in shape for dark matter and radial bins, but 
the amplitude in this case is 0.1. 

If we bin the plane n — a with IMpc/h, the same theory 
error form in Ea. (|A4|l still works on all scales. So does the 
JK, at low a. But in general errors are higher than when us- 
ing a wider bin (smaller bin means smaller number of pairs). 
This increase is compensated by a decrease in the covariance, 
which becomes practically zero for all elements outside the 
diagonal. 
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Figure A6. Left: Monopole ^(s) with JK errors for MICE simulations (black, lower line) and for real LRG (red, upper line) using a bin 
of 5Mpc/h. Right: Diagonal error for MICE dark matter simulations (JK with dispersion solid black line, MC dashed black line) and for 
JK errors in real LRG data (red for Nr = 10 and blue for A'^^ = 20) 



A3 Errors in the monopole ^(s) 



The error in the monopole is not easy to predict theoret- 
ically, but in this case JK error works very well for both 
variance and covariance on all scales. In the left panel of 
Fig |A6l we plot the monopole when we bin the data with a 
separation 5Mpc/h: black for the dark matter mocks and 
red for LRG data (which will be presented in more detail in 
section |4]|. In the right panel, we see the difference between 
the mean JK error and its dispersion (solid line with errors) 
and the MC error (dashed line), and over-plotted the JK 
error for the real LRG (color). As can be seen, the error for 
LRG is larger than the one in simulations, since the signal 
is also higher, due to an overall bias b in the amplitude. The 
small difference between blue and red line in Fig |A6l is due 
to the number of the particles in the random catalog (A_r), 
which does not change the estimation of the error in this 
case. Results are very similar for group mocks, with higher 
amplitudes, as expected. The important point here is that 
the JK error works well at all scales. 



The normalized covariance, in Ea. (|A3|l . for the monopole 
is plotted in Fig |A7l The result is not completely smooth for 
the MC case (top panel), because we need more than 216 
simulations to have a smooth result. But we see that it has 
the same shape than JK covariance, which we will use to fit 
our LRG data. JK covariance is smoother than the MC be- 
cause we have taken the mean over 216 x 63 realizations (216 
mocks X 63 JK zones), which seems enough to converge. The 
bottom panel shows JK errors in the real LRG data, which 
is quite similar to the mocks, but looks noisier because this 
is just a single realization. 

We have also binned the monopole ^(s) with IMpc/h. 
JK error also works in this case, and we also see that the 
error is higher here than in the 5Mpc/h bin, while the co- 
variance is smaller, practically equal to zero. 




Figure A7. Normalized covariance of MC (top) and JK (mid- 
dle) for the monopole ^(s) in MICE group mocks (for binning of 
5Mpc/h). Bottom panel correspond to JK error in the real LRG 
data 
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Figure A8. Top Left: Q(s) with errors for MICE dark matter mocks (black higher signal) and for LRG (red lower signal) using a bin 
of 5Mpc/h. Top Right: Diagonal error for MICE simulations (JK with dispersion solid black line, MC dashed black line) and for LRG 
data (red and blue for A^^ = 10 and A''^ = 20). Bottom: Same as top with a bin of IMpc/h. 



A4 Error in the quadrupole Q(s) 

For the quadrupole, the JK error also works quite well, as 
we can see in Fig lASI for 5Mpc/h (top panels) and IMpc/h 
binning (bottom panels). Again here , the solid line with 
errors shows the JK error and its dispersion, the dotted line 
is the MC error, and the color and lower line is the JK error 
in LRG data. The error for LRG data, in the quadrupole, is 
lower than the simulations one. This is because the error here 
is proportional to the signal, and in this case, LRG signal is 
lower since Q(s) depends on /3 which depends inversely on 
bias b. Results for group mocks are in good agreement with 
data and yield similar results for the comparison on JK and 
MC errors. 



A5 Validity of the models 

We also use the simulations to test the methods that we will 
apply to real data (LRG) in the following sections. We want 
to study if we can recover the parameters that were input 
in the simulations. 



A5.1 Recovering the real-space correlation function 

We first estimate the 2-point real-space ^(r) and redshift- 
space ^(s) correlation functions and compare ^(r) to the 



input model that describes the simulation (see top panel in 
Fig |A9|) . We find very good agreement. We look at the ratio 
5(s)/^(r) between the redshift space and real space, which 
should be a function of the distortion parameter /3 at large 
scales (and for the distant observer approximation), as in 
Eq. (|23|l . In bottom panel of Fie lA9l we have plotted the mean 
ratio over the simulated mocks, with its error, and over- 
plotted in red the expected value for the ratio at /3 = 0.62 as 
in the input model. We also plot the mean(5(s))/mean(^(r)) 
in dashed line (a biased estimator), which is similar to the 
mean(^(s)/^(r)) until 30Mpc/h. It seems to converge below 
lOMpc/h and starts to fail around 30Mpc/h, where errors 
rapidily increase, which is in agreement with other analysis 
(see Fig. 13 in Hawkins et al 2003). We could obtain f3 from 
the ratio ^(s)/^(r), but it is difficult in real data since we 
do not have direct information of the real-space ^(r), only 
through integration of the anisotropic ^{a, n) through the 
line-of-sight (see i^2.3|l . The advantage of using this ratio is 
that it converges to a constant value at small scales. But we 
can only use this below 30 Mpc/h with real data because the 
recovered real-space correlation function ^(r) starts to fail 
on larger scales (see Fig[l|. We have fitted these scales with 
the simulations and the results for /3 are the ones expected. 
However, we have checked that it is better to obtain /3 from 
the quadrupole Q(s), since it is more independent of other 
cosmological parameters and can be trusted to larger scales. 
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Figure A9. Top panel: Real-space correlation function ^(r) 
(black dots) with its model (solid line) and redshift space corre- 
lation function ^(s) (red dots), which is f (r) biased by a constant 
Kaiser value for large scales. Bottom panel: Mean i{s)/^{r) over 
MICE dark matter mocks with errors (on the mean, scaled as 
l/Vnsim), and Kaiser prediction for large scales for /3 = 0.62, 
corresponding to the input model (red). In dashed, we have plot- 
ted mean(5(s))/mean(^(r)), a biased estimator, which is different 
than the solid line mean(^(s)/^(r)). The ratio is constant from 
8Mpc/h up to 30Mpc/h, where errors start to blow up. 
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Figure AlO. Dispersion in the pairwise velocity distribution CTi, 
in the MICE simulation as described in Eo. 11511 when we change 
the distance between particles 




s(Mpc/h) 

Figure All. Top panel: Mean best fit j3 - av to MICE 
quadrupolc Q(s). The red dot shows the value for the input f} and 
cr„ at large scales. Contours are Ax^ = 1,4 and 9 using MC errors 
for a single mock. Bottom panel: Q(s) with MC errors (points), 
and best model (red). We also plot in blue the quadrupole for 
one of the mocks, in order to see the scatter from bin to bin. We 
recover the input /3 = 0.62, but Oy is slightly biased to higher 
values, probably because we are obtaining an effective cr^ which 
also accounts for the values at lower scales. 



A5.2 Fitting of the quadrupole Q(s) to obtain (5 and o-u 

We can calculate the multipoles of ^(tt, a) to decompose the 
anisotropy between the LOS and the perpendicular direc- 
tion. We have tested with models that the monopole ^(s) 
and quadrupole ^(s) and even the combination £,{s)/^{s) 
depend strongly not only on (3 and o^ , but also on other pa- 
rameters like the shape of the correlation function (ie flm, 
fit. and Us), the non- linear bias and the overall amplitude. 
On the contrary, the reduced quadrupole Q(s), defined in 
Ea. (|17p only depends strongly on (j„ and (3, but not on the 
bias. So when using the quadrupole we do not need an ex- 
pression for the non-linear bias or a prior on flm (the shape 
of ^(r)) to extract the /3 information. We can fix these pa- 
rameters, and only change /3 and o-^. The asymptotic value 
of Q(s) for large scales in Ea. (|18p is a function of /3, but 
it is always biased to higher values because of the random 
velocities. So it is dangerous to obtain (3 from this asymp- 
totic approximation. We instead generate our model for Q(s) 
based on the multipoles of ^{a, n) in Eq. (fT6)) and Eq. (fT7)) . 
For the dark matter mocks we have that /3 = 0.62 (see 
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Figure A12. Correlation function f (tt, <t) in different simulations: MICE dark matter mocks with 6 = 1 and fS = 0.63 (left); group mocks 
with b = 1.9 and /3 = 0.25 (middle); and group mocks with b = 3.0 and /3 = 0.23 (right). The contour colors are from the simulations. 
We have overplotted in solid lines the corresponding levels for models with the same input parameters as the simulation in each case. 



^XJ and av is approximately 400km/s on large scales, but 
higher when approaching scales smaller than 10 Mpc/h (see 
Fie: IA10( ). We fit the quadrupole obtained for each mock sep- 
arately. The mean over all the mocks gives the correct exact 
value of /3 = 0.62 and a value of (t„ ~ 450Km/s which is 
slightly larger than the asymptotic large scale value of 400 
Km/s. This is probably because we are obtaining an effective 
a-u which also accounts for the values at lower scales which 
are larger, as shown in Fig lAlOl But at 1 — cr (2 dof) a^ ~ 
450Km/s is consistent with the asymptotic value at large 
scales, which is plotted as a red dot in Fig lAllI In Paper-II 
of this series we show how to find the exact value for (t„ as 
a function of scale. 

In the top panel in Fig lAllI we show the mean contour 
P ~ (Jv which is the average over the individual contours in 
each simulation, obtained from Q{s). In the bottom panel, 
we show the best fit model over-plotted over the mean Q(s) 
with the MC errors. The MC errors here correspond to a 
single mock, while errors in the mean value are viV times 
smaller (N=216 in our simulations). In the bottom panel of 
Fig lAllI we plot the errors corresponding to a single mock. 
We find similar results for group mocks, which have differ- 
ent values of 13. We conclude that the model used for the 
quadrupole, assuming a constant Ov, is perfect to obtain 13 
and an effective av which does not have to correspond nec- 
essarily to the large scale value, but rather to a combination 
of (T„ at different scales. 



A5.3 Fitting of large scales: Q.m and amplitude 

From large scales, we will obtain the shape of the linear 
correlation function ^(Tr, ct). Here we will only fit fim and 
the amplitude of the correlation. We have fixed Us = 0.98, 
Q.b = 0.045 and h = 0.72 following recent results of WMAP, 
SNIa and previous LSS analysis Q 



see |http://lambda.gsfc.nasa.gov/product/map/dr3/parameters.cfm] 



As we have seen in i|2.3l we can not rely on the projected 
correlation function at large scales, or the recovered real 
space correlation function, so the way to extract information 
about the shape of the real space correlation function is 
directly from ^(cr, vr). We use models that vary with Q,m, 
linear amplitude, a^ and (3. The linear amplitude Amp refers 
to the factor b{z)ag,, where b(z) is the bias at redshift z. 
b{z) and as are completely degenerated in the correlation 
function, also with the growth factor D{z). In fact, what we 
obtain from observations is b{z)af!,D{z), but D(z) is known 
for each cosmology and the median redshift of the slice. We 
have checked with simulations that we can use linear bias 
in the scales that we are interested. We also use the priors 
(T„ and (3 found previously from the normalized quadrupole, 
although this does not make a large difference. Then, we 
need to fit Jim and Amp. The best place to fit them, so 
that is independent on scale, is for intermediate scales, from 
20Mpc/h to 60Mpc/h, to be safe from non-linear bias. We 
also stay away from the line-of-sight (LOS) where fingers of 
God can alter the information, cutting an angle of 30-40 deg 
away from the LOS when doing the fit. We marginalize over 
(3 — av and obtain a contour for fJm — Amp which is in very 
concordance with the input values. Although a^ obtained 
from Q(s) is just an effective value, simulations show that 
we recover the correct value of f2m- We can understand it, 
because av does not play an important role in our selected 
range of scales in ^(a, tt). If we extend to larger distances, 
the value for Q.m is biased slightly to lower values, probably 
because of the wide angle eflect. This is important for an 
angle of 10 deg between galaxies (Matsubara 2000), which 
is more than 70Mpc/h for our worst case (z=0.15), in the 
perpendicular direction a. 

A5.4 Model for the 2-point correlation function ^{a, n) 



Finally, in Fig |A12l we plot ^{a, tt) at large scales for the 
mean over MICE mocks (colors) . We over-plot in solid lines 
the best model (using Ea. (|14|l ). In all cases, we use an effec- 
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tive ffv = 400km/ s. Note that we use square bins (pixels) of 
5Mpc/h side at all scales and this dilutes the FOG at small 
a. The model works very well compared to simulations. We 
can observe in the figures that the overall amplitude is dif- 
ferent for each simulation, because it is proportional to the 
bias. The same bias modifies the distortion parameter /3 (in- 
verselly proportional to the bias) which changes the shape 
of (,{n,a). The larger the value of l3 the greater squashing 
effect in the radial (tt) direction. Also note the closed con- 
tours in the radial direction between tt = 50—100 Mpc/h (in 
dark blue) . They correspond to a region of negative ampli- 
tude caused by the squashing (Kaiser) effect. The larger the 
value of /3 the larger the region with negative values. This 
region is sorounded by the BAO ring at about 100 Mpc/h. 
It is remarkable how well this is followed by simulations, 
which indicates that this is a very significant feature. We 
have shown before that data also follows this feature. 

Simulations also show a good agreement with the BAO 
peak, including some detection in the radial direction which 
is enhanced by noise. This will be studied in more detail in 
Paper IV of this series. 

We see that the obtained correlation (colored) differs 
slightly from the distant observer approximation theory (lines) 
at large a and tt ~ 0. The redshift space correlation distor- 
tion in real surveys, which are not located at infinity, de- 
pends on n and a, but also on the angle between galaxies 9 
and the angle 7z between the direction LOS (at 9/2) and the 
vector which goes from galaxy 1 to galaxy 2 (following the 
notation used in Matsubara 2000). In the distant observer 
approximation we assign to the angle between galaxies the 
value ^ = 0. Matsubara (2000) has studied the differences 
between the real correlation and the approximation for dis- 
tant observers. In general, the approximation is good for 
angles 8 below 10 deg, which include all the zone we are us- 
ing for our analysis in LRG data (the worst case comes from 
the closest galaxies at z=0.15, where a ~ 80Mpc/h corre- 
sponds to 6' = lOdeg). The correlation also depends on the 
7z angle. We can see in Fig. 9 of Matsubara (2000) a com- 
parison between the real correlation, which depends on the 
distance between galaxies and both angles described above, 
and the distant observer correlation function, used in this 
work, which depends only on n and a. Each position in the 
2 dimensional tt — ct is a mixing of different 9 and 7z, because 
there is a range in redshift, but we can explain qualitatively 
the lack of power of the observed correlation respect to the 
distant observer approximation theory at large a and tt ~ 0. 

We have also calculated the correlation function limit- 
ing the angle between galaxies to see in practice how ^(tt, a) 
is affected by wide angle effects (see Appendix B5). As we 
increase the restriction, we see how the a direction recovers 
power in better agreement with the model predictions. In 
real data, we rather use all pairs given that this is quite a 
small distortion. 
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Figure Bl. Selection function for data (black histogram) and 
smoothed for the random catalogs with a bin in redshift of z=0.02 
(red), z=0.05 (green) and z=0.08 (blue) 



APPENDIX B: SYSTEMATIC EFFECTS 

In this section we look for possible systematic errors that 
could be imprinted by the radial mask through the line-of- 
sight, the angular mask, or the selection of LRGs . These sys- 
tematic effects are typically more important on the largest 
scales, where the correlation function becomes smaller. In 
Fig llQI we show that for some slices there is extra power in 
f(s) at scales s > 130Mpc/h, something which is not ex- 
pected in the models. This is particularly evident for the 
slice in redshift z=0.3-0.4, where the signal-to-noise is small 
(Amp and bias 6 are smaller in Table [Sjl and we do not have 
a significant peak detection. This extra power could be due 
to sampling variance but could also be caused by some sys- 
tematic errors in the data or the way we analyze it. The 
extra power seems to be important at scales larger than the 
baryonic peak, but we test here if it could also have some 
effect over the peak location. This extra power have also 
been detected in other analysis of SDSS LRG around simi- 
lar scales (Blake et al 2007; Padmanabhan et al 2007). We 
will focus here on systematics related to the way we have an- 
alyzed the data, rather than systematics in the data itself. 
The later has been explored in detail elsewhere (Eisenstein 
et al 2005). 

Bl Radial Selection 

First, we test the radial selection function that we use for 
the random catalogs. If we use exactly the same radial den- 
sity distribution A'^(2;) as selection function of the data, we 
suppress the radial modes in the tt direction. In Fig lBll we 
can see the differences between different smoothing windows 
in the data selection function, and in Fig |B2l and Fig |B3l the 
redshift-space correlation function for these three smoothing 
bins. 

We do not see any significant difference between the 
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Figure B2. Redshift-space correlation using different redsliift 
smoothing in the random selection function of 0.02 (red), 0.05 
(green) and 0.08 (blue) 



three cases indicating that our analysis is robust with respect 
to the radial selection. 



B2 Weighting 



In Fig |B4l we compare our results in the monopole (red with 
shaded region as errors) to that of Eisenstein et al (2005) (in 
black line with errorbars) for all the sample. Our result is 
consistent with this previous estimate despite the increase in 
the DR6 area and the difference in the selection. However, 
at larger scales than the baryonic peak, we observe some 
extra-power in our estimation. We wonder if the difference 
is just due to sampling, selection or to the way we estimate 
the correlation function (we do not include the same weight- 
ing as Eisenstein et al 2005) 

We have also calculated the correlation function using 
a weighting scheme, as the one explained in Eisenstein et 
al (2005). We weight the sample using a scale-independent 
weighting that depends on redshift. When computing the 
correlation function, each galaxy and random point is weighted 
by 1/(1 + n(z)Pw) (Feldman et al 1994) where n(z) is 
the comoving rmmber density and Pw = 40, 000/i~"^Mpc^. 
We do not allow Pw to change with scale so as to avoid 
scale-dependent changes in the effective bias caused by dif- 
ferential changes in the sample redshift. We choose Pw at 
k ~ O.O^hMpc'^ as in Eisenstein et al (2005). At z < 0.36, 
nPw is about 4, while nPw ~ 1 at z = 0.47. 

In Fig lBSI we can see the comparison between the corre- 
lation function estimated without weighting (solid line) and 
with weighting (dashed line). Contrary of what we were look- 
ing for, the extra power is higher in the weighting scheme, 
which makes sense, since we are now giving more impor- 
tance to the higher redshift pairs, which have a larger bias 
(see below). 
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Figure B3. ^(Tr, a) correlation using different redshift smoothing 
in the random selection function of 0.02 (top), 0.05 (middle) and 
0.08 (bottom). 
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Figure B4. Estimation of tlie redshift-space correlation function 
for LRG. In black line with JK errors correspond to the result by 
Eisenstein et al. 2005. Over-plotted in red is our result, where the 
shaded region corresponds to our error estimate. 
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Figure B5. Our previous estimation of the redshift-space cor- 
relation function for LRG (solid) compared to estimation using a 
weighting as explained in the text (dashed) 
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Figure B6. Our previous estimation of the redshift-space corre- 
lation function 5(7r, c) for LRG (top panel) compared to estima- 
tion using a weighting as explained in the text (bottom panel) 



In Fig lBGI we have plotted the anisotropic correlation 
function ^(71,0-) for both cases, showing only slight differ- 
ences. It is apparent that the monopole is more sensitive that 
^(71,0-) to displace systematic differences. We will therefore 
concentrate on results for the monopole from now on. 

B3 Angular Selection 

We next look at the angular selection function, the mask. 
First, we construct a different mask than the original by 
using a Healpix map (Gorski et al 2005), with nside=64 
(with pixels of area ~ 0.8 sq dog). In Fig |B7l we plot the 
distribution function for the number of galaxies per pixel. 
The pixel is large enough to distinguish between real empty 
pixels and artificial ones. We only include pixels that have 
more than n galaxies, where n=2,6 and 10. As we can see 



in the distribution plot, if we include more than 2 galaxies, 
we are probably including artificially void zones, which will 
create extra power and pencil beams at the direction line-of- 
sight, while when we include only the pixels with more than 
10 galaxies, we are probably excluding some real voids. In 
this second case, the density is higher than the real one, 
so the density contrast is lower. We can see these effects in 
Fig lBSI for the correlation function in redshift space. We see 
that the correlation is lower when we increase the minimum 
number of galaxies (from red to blue), but we can always 
see the baryonic peak and the three lines have the same 
shape approximately. Our results are similar to the case n 
> 6, which does not include artificial voids and does not 
eliminate real voids. 

The spectroscopic survey of SDSS is observed using cir- 
cular plates, which contain about 600 fibers each to take 
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Figure B7. Distribution function of number of galaxies per pixel 
(nside=64 using the package Healpix) 



Figure BIO. Redshift-space correlation function for our mask 
(black), maskplatel as indicated in the text (red), maskplate2 
(dotted red) and north stripe (blue) 
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Figure B8. We now define our catalog by including all the pixels 
that have more than n galaxies in the Fig |B7l (n>2 red, n>6 green, 
n>10 blue) and plot the redshift-space correlation function 



spectra. Targets are selected from the photometric survey, 
although the spectroscopic survey is not exactly the same as 
the photometric one. There are some plates that have not 
been observed properly due to known problems, which are 
explained in the web (http://www.sdss.org). We have ex- 
tracted from our previously calculated mask all the galaxies 
that are not laying inside "good" plates (maskplatel). In 
Fig |B9l we can see the plates (black circles) and the galaxies 
(red dots). Moreover, we can also eliminate all the galaxies 
that lay inside a bad plate to ensure that we are taking only 
the really good ones (maskplate2). This second mask reduces 
the number of galaxies significantly, and the correlation is 
then a lot noisier, but we show here the results. In Fig lBlOl 
we have plotted the redshift-space averaged correlation func- 
tion for the new mask based on spectroscopic plates. Results 
are very similar to our previous result, in black. Moreover, 
we over-plot the result for the north stripe of SDSS, which 



contains the most significant part of the survey, in blue. The 
anisotropic redshift-space correlation function is really simi- 
lar in all these cases, so that we do not find it useful to plot it. 

We have also tried to extract all the plates that have a 
large number of galaxies, which have big superclusters, since 
it is known that big clusters can bias the correlation func- 
tion for regular galaxies (Baugh et al 2004). The result is 
not significantly modified which could have been expected 
since LRG in SDSS DR6 cover a much larger volume that 
2DFGRS or SDSS for regular galaxies. 



B4 Volume limited samples 

In volume limited surveys, we can estimate the correlation 
function with a pixelization scheme. The correlation func- 
tion can then be estimated as: 



C(s) = ^ SG{Si)SG{Sj)/Npairs{s) 



(Bl) 



where Sq = No/ < Nq > — 1 and the sum extends 
to all pairs i,j separated by a distance s ± As. We have 
taken a volume limited part of the selected galaxies, with 
redshift z=[0.15,0.38] and absolute magnitude Mr=[-22.5,- 
21.5] and have calculated the correlation using this method. 
This is the same sample used in Paper IV of this series. We 
have also calculated a similar selection by using the tradi- 
tional method with a random catalog. Results are plotted 
in Fig lBllI dotted for the pixel method and solid for the 
randoms. This is a good test to validate our results, since 
both methods are quite different. As we can see in the fig- 
ure, the two estimations are very similar. The slight shift in 
amplitude at small scales is due to the pixel smoothing, a 
cube of 10 Mpc/h on the side for this case. This shift can 
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Figure B9. SDSS DR6 survey with our selection of LRG galaxies as red points and the spectroscopic "good" plates as black circles 



Table Bl. Slices in the plane Mr 
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Figure Bll. Redshift-space correlation function for a volume 
limited slice selection of LRG galaxies with z=[0.15,0.38] and 
Mr=[-22.5,-21.5]. Solid line shows the result when using a random 
catalog. Dotted line shows a new method based on pixelization 
which validated the previous results 



be modeled with a top-hat window function reproducing the 
exact amplitude difference. There is a very good match using 
smaller pixel size. Also note how these estimations, based on 
a volume limited sample with about half of the LRG galax- 
ies, agree quite well with previous results (ie in Fig lBlOJ I 
which includes all galaxies at a price of a more complicated 
selection function. Despite all this difference, we see no sig- 
nificant change in ^(s) at the largest scales. 

We have also divided the catalog in different volume 
limited slices as indicated in table IB4I plotted in Fig |B12l 
The most and least luminous slices (in the bottom right and 
top left) contain fewer galaxies and have too much noise at 
the baryonic peak, but we see that the red and pink slices, 
at intermediate redshifts and magnitudes, contribute to the 
extra-power that we see at large scales (Fig |B13l with the 
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same colors as Fig |B12|l . Note how these subsamples cover 
a sharp feature a z=0.35 in the selection function (eg see 
Fig lBlf) which is clearly an artifact of the selection function 
because of the double selection cut in LRG. This feature 
is diluted when we consider volume limited samples as we 
discard the less luminous galaxies Mr > —21.5 which con- 
tribute most to this peak. As shown above, the correlation 
5(s) for a volume limited sample in the range Mr=[-22.5,- 
21.5], shown in Fig lBlll is in excellent agreement with the 
one for the whole sample in Fig |B4l Thus we conclude that 
our measurements are quite robust and are not affected by 
the selection cuts. 

We have not plotted larger scales for the other slices be- 
cause they are quite noisy. However, at intermediate scales, 
we can see clearly the bias due to different intrinsic luminos- 
ity, more biased when more luminous, although bias seems 
to be independent on scale, in the scales used for our anal- 
ysis. 



B5 Wide angle approximation 

We have calculated the correlation function limiting the an- 
gle between galaxies to see if wide angle effects that are ap- 
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Figure B12. We have divided tlie catalog in different approx- 
imately volume limited slices as indicated in table IB4I Here we 
over-plot the slices in the plane Mr — z 
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Figure B13. Redshift-space correlation function for the different 
slices plotted in Fig |B12l 
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Figure B14. Redshift-space correlation function ^(Tr, a) includ- 
ing all the galaxy pairs (top panel) and limiting the angle between 
galaxies to S < lOdeg (bottom panel) 



parent in ^(tt, a) disappear. In Fig |B14l we see the anisotropic 
^(tt, a) without limits in the angle (top panel) compare to 
the one only accepting galaxies with 9 < lOdeg (bottom 
panel). As we increase the restriction, we see how the a di- 
rection recovers power, which explains the lack of power we 
saw in the a direction when angles are too big to apply the 
distant observer approximation. The angle between galaxies 
explain part of the distortions due to wide angle effects, spe- 
cially the ones that are concentrated at small vr and large 
a, which affect the first slice considered, from z=0. 15-0.3. 
The angle 7z, between the direction LOS (at 9/2) and the 
vector which goes from galaxy 1 to galaxy 2 (following the 
notation used in Matsubara 2000), is also important and can 
also imprint some modifications at larger n and large a. 



probably the result of sampling fluctuations and shot noise. 
For redshifts around z=0.35 we find extra noise on large 
scales, probably due to an artifact in the selection function. 
In any case, none of the systematics we have explored modify 
the peak detection in a significant way. 



We conclude that the measurements on the BAO scale 
are quite robust and the extra power at the largest scales is 



© GOOD RAS, MNRAS 000, 000-000 



